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ABSTRACT 

The physical characteristics of Cepheus E (Cep E) 'embedded' outflow are 
analyzed using ISOCAM images in the v=0-0 S(5) 6.91 //m and S(3) 9.66 
//m molecular hydrogen lines. We find that the morphology of the Cep E outfiow 
in the ground vibrational H2 lines is similar to that of the near infrared v=l-0 2.12 
//m line. At these mid-IR wavelengths, we do not detect the second H2 outflow 
which is almost perpendicular to Cep E 2.121 fj,m flow or traces of H2 emission 
along the second ^^CO J = 2 — 1 outflow at ~ 52° angle, down to a surface 
brightness of 12 - 46 /xJy/arcsec^. 

We do detect at 6.91 fj,m. the likely source of the main H2 and CO outflows, 
IRAS 23011+6126, and show that the source is easily seen in all IRAS bands using 
HiRes images. The source is not detected at 9.66 /xm, but we think this agrees 
with the interstellar extinction curve which has a minimum at ~ 7 /im, but rises 
a ~ 9.7 fj,m due to the strong absorption silicate feature, enhanced in this case by 
a cocoon surrounding the Class object. This idea is supported by our models of 
the spectral energy distribution (SED) of the central object. The models assume 
that the main source of opacity is due to bare silicates and our best flt for the SED 
yields a total mass of envelope of 17 Mq and a dust temperature of 18 K. 

Subject headings: stars: formation — stars: pre-main-sequence — ISM: jets and 
outflows — infrared: ISM: lines and bands 
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1. Introduction 

The Cepheus E (Cep E) outflow is an excellent object to study the relationship between 
the physical properties of the optical stellar jets and those of embedded outflows. The 
existence of a molecular CO flow originating in that region was flrst indicated by Fukui (1989) 
in his catalog of molecular outflows, but it was the K' image of Hodapp (1994), from his 
imaging survey of molecular outflows, that brought Cep E to the forefront. The south lobe of 
Cep E is observed optically in Ha and [SIX] AA6717/31 (Noriega-Crespo 1997), and it has been 
named Herbig-Haro (HH) object 337 (Devine, Reipurth & Bally 1997). The spectra of this 
lobe displays one of the lowest ionizations measured in HH objects, with a ratio [SII]/Hq;= 
8.8 ±0.2 (Ayala et al. 1998). 

The Cep E is bright in H2 emission lines at 2/xm and is likely to be driven by the IRAS 
23011+6126 source, presumably a class Class object (Eisloffel et al. 1996). Because of its 
complex morphology and strong H2 emission, Cep E has been used to gauge the reliability 
of 3-D molecular jet models (Sutter et al. 1997). There is also clear evidence for a second 
outflow in H2 almost perpendicular to the main flow emission, as indicated by some faint H2 
knots (Eisloffel et al. 1996). Perhaps more surprising is the presence of a second molecular CO 
outflow, detected in both ^^CO J = 2 — 1 and ^^CO J = 1 — molecular lines (Ladd & Hodapp 
1997; Ladd & Howe 1997), at an angle of ~ 52° from the main CO flow, i.e. apparently 
unrelated to the second H2 outflow. The ^^CO emission indicates terminal velocities of 80 
km s~^ and —125 km s~^ in its north and south lobes respectively, which implies a dynamical 
age of ~ 3 X 10^ years, at a distance of 730 pc (Eisloffel et al. 1996). The absence of a second 
source using intcrfcromctric observations at 2.65 mm suggests that both outflows arise very 
close to the position of IRAS 23011+6126 source (Ladd & Howe 1997). 

In this study we present new ISOCAM images of the Cep E outflow taken in the ground 
vibrational level (v=0-0) H2 lines S(5)at 6.91 fj,m and S(3) at 9.66 fj,m. The motivation of 
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this work was to use the ground vibrational H2 hnes, which are predicted to be much brighter 
than the v=l-0 S(l) Une (Wolfire & Konigl 1991), to study the excitation across the outflow 
(using the S(3)/S(5) ratio), and to try to detect the second H2 outflow and the central source 
at mid-infrared wavelengths, as well as to search for traces of H2 emission along the second 
CO outflow. These observations are complemented by HiRes IRAS images. 

2. Observations 

The Cep E outflow has a projected size onto the plane of the sky of ~ 1', so the ISOCAM 
images were obtained using a 2x3 CVF raster map with a 6" FOV pixel scale and 30" steps. 
Four CVF fllters were used; two very close centered on the v=0-0 S(3) 9.665 /im and S(5) 
6.909 B.2 lines respectively (see below); plus two nearby continuum CVF steps at 9.535 
/im and 6.855 yum. We selected the S(3) 9.665 //m line, despite the fact that its wavelength is 
right in the middle of the strong silicate absorption feature at ~ 9.7 fj,m, because the line is 
expected to be strong and should help to constrain the depth of the silicate absorption in the 
models of the spectral energy distribution (SED). 

The ISOCAM data was deglitched using the Multi-resolution Median Transform Method 
(CIA). The detectors transient effects were treated using the IPAC Model (Ganga 1997). The 
target dedicated time (TDT) was of 2728 seconds, with twelve stabilization time steps on hne 
and ten in the continuum prior to the on target observations. The TDT was spent with 2/3 of 
the time on line and 1/3 in the continuum, this means approximately 900 sees in each of the 
H2 lines. The ADU fluxes were transformed into calibrated fluxes using the upgraded values 
of the system response (e. g. ISOCAM Observer's Manual, Tables 12-17). These values are 
given in ADU/scc/mJy/pixcl and correspond to 125.52 (step 227 9.660 /xm), 65.80 (step 331 
9.535 /xm), 122.12 (step 21 6.911 /xm) and 121.686 (step 22 6.855 fxm). The FWHM of four 
the fllters are (9.660, 9.535, 6.911, 6.855) fj,m = (0.27, 0.22, 0.17, 0.17) pim; and we used an 
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integration time step of 2 sees and a ADC gain of 2. 

We present also for comparison a near infrared image at the v=l-0 S(l) 2.12 iim obtained 
at the 3.5m Apache Observatory with a 256x256 array at f/5 with a 0.482"per pixel scale. A 
complete analysis of the imaging and spectroscopic near infrared data is presented elsewhere 
(Ayala et al. 1998). HiRes IRAS images are presented also to support these observations. 
These images have one degree field of view with 15" pixels and can reach a spatial resolution of 
~ 1'. The HiRes images have been processed using Yu Cao algorithm (see e. g. Noriega-Crespo 
et al. 1997). 



3. Results 
3.1. Morphology 

The grayscale images of H2 at 6.91 ^m and 9.66 /um are presented in Figures 1 and 2 
respectively. From these images is evident that the morphology of these different molecular 
hydrogen lines is very similar. It is also clear that the main difference between them is the 
central intensity peak at 6.96 /im. This intensity peak coincides with the position of the IRAS 
23011+6126 source, as determined using interferometric observations by Eisloffel et al. (1996), 
i. e. a(2000) = 23h 03m 13.0s, 5(2000) = 61° 42' 26.5". The IRAS 23011+6126 source is also 
detected in the continuum frame at 6.855 //m, indicating that the emission is not dominated 
by the excited H2 emission. The Figure 3 shows the 6.91/im image with an overlay of the H2 
2.12 //m emission, and once again the morphology of both molecular lines is very similar. 

One way to understand these observations is in terms of the interstellar extinction law 
(see e. g. Mathis 1991), which it has a minimum at ~ 7 ^m and the rises to a peak at ~ 10 
/im, where this maximum is due to the silicate absorption feature. If the IRAS 23011+6126 
is embedded in a dusty envelope, as expected for a Class source, then the absorption by 
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silicates will be even larger and enough to swamp the H2 emission at 9.66 //m. The IRAS 
23011+6126 source appears in all the IRAS bands, as is shown in Figure 4, which displays 
IRAS HiRes maps at 12, 25, 60 and 100 /^m centered on the Cep E source. 

We recall that another two outflows have been detected around the Cep E source. One 
outflow is observed in the 2.12 /im H2 emission and is almost perpendicular to the main H2 
flow. The second one is detected in the ^^CO J = 2 — 1 transition and is centered in the 
IRAS 23011+6126 source (±2.3"), with an orientation of ~ 52° with respect the main H2 flow 
and a scale of ~ 4' (Ladd & Hodapp 1997). We don't detect in the H2 lines at 9.66 /xm or 
6.91 fj,m any signatures of the second H2 outflow nor of the second CO outflow observed at 
millimeter wavelengths. For the faint 2.12 //m H2 outflow this perhaps is not surprising since 
the S(3) 6.91 iim line would have needed to be ~ 70 times stronger than S(l) 2.12 /xm line 
(based on simple C-type shock models (e.g.Smith 1995)) to overcome its very low surface 
brightness in comparison with the brightest regions. We estimate that the RMS noise levels 
of our ISOCAM images, based on measurements of the background are approximately 46 
//Jy/arcsec^ at 9.66 /xm, and 12 //Jy/arcsec^ at 6.91 //m. The values of the minimum contours 
in Figures lb and 2b (continuum subtracted) are 0.4 and 0.2 mJy/arcsec^ for the 9.66 fj,m and 
6.91 yum lines respectively, so a factor 70 times lower in brightness for the faint H2 outflow is 
at the noise level. There are not traces in the IRAS HiRes images of the faint H2 flow nor 
of the CO outflow. The HiRes images do not resolve either the second source near IRAS 
23011+6126, as is shown in Figure 4. The 12 and 25 //m images in Figure 4 display only point 
sources, and although there is some faint emission at 60 /xm and 100 /xm at a PA~ —45°, this 
is probably be due to the diffuse emission from the nearby sources. The existence of a second 
source at 2' SW from the IRAS 23011+6126 source has recently been conflrmed by L. Testi 
(1998) with OVRO observations at 1.3 mm. 
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3.2. Fluxes 

As mentioned before the v — — S(3) and S(5) H2 were selected for our observations 
because plane parallel molecular shock models predict, for the conditions encountered in HH 
objects, that these lines could be ~ 10 — 100 stronger than the v — 1 — S(l) 2.12 /im line if 
they are coUisionally excited. Molecular shock models specifically calculated for HH objects 
(Wolfire & Konigl 1991) considered as typical parameters e. g. a shock velocity of Vg =25 
km s~^, an initial preshock gas density of uq = 10^ cm~^ and a preshock magnetic field of 
Bq — 30 yuG, a molecular hydrogen abundance ~ 0.5, an atomic hydrogen abundance 
^ ~ 3 X 10~^, a neutral gas temperature of T ~ 2000 K and an electron temperature 
Te ~ 3000 K. The point of ennumerating some of these parameters (there are a few more) is 
to illustrate the difficulty in comparing shock models with observations and to stress that the 
models should be taken as a guide, since the physical conditions across a shock front change 
and are more complex. 

For the input parameters of the shock models mentioned above, the expectation is 
that the ratios of the H2 lines should be '^^^ ~ 157 and ^Ip g^^j ~ 28, or essentially 

=5.6. Figure 5 shows the ratio of our S(5) to S(3) H2 images, which indicates the the 
ratio across the outflow is nearly unity. The ratio is constant except for the region around the 
IRAS 23011+6126 source, which is not detected in the S(5) 9.66 /xm image, and at the edge 
of the outflow, which is probably due to an artifact produced by the mismatch between the 
elliptical Gaussian function use to smear out the S(3) 6.91 /im image, and the true shape of 
the first Airy ring produced by the ISOCAM optics. 



Nearly constant distributions across the outflow have been also measured for the ^"^ g^^^ 
and ratios (Eisloffel et al. 1996). This behavior of the ratios is very difficult to 

explain with simple shock models or shock geometries, and in the case of the v=2-l S(l) and 
v=3-2 S(3) lines ratios, the best results were obtained with C-type bow shocks. Such C-type 
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bow shock model needed shock velocities of ~ 200 km s ^ and preshock densities of ~ 10^ 
cm~^ (Eisloffel et al. 1996), which are larger than what is measured in most optical outflows. 

A constant v=0-0 S(5) to S(3) ratio also indicates the lack of a steep extinction gradient 
between the north and south lobes, and that the extinction is mostly important around the 
IRAS 23011+6126 source. This is interesting because we know that the south lobe is visible 
at optical wavelengths (Noriega-Crespo 1997), while the north lobe is not. 

Finally, we have measured the flux of the IRAS 23011+6126 source at 6.855 iim , set an 
upper limit for the 9.535 fj,m flux (which may be still affected by the silicate feature) and 
measured the IRAS fluxes from the HiRes images. In Figure 6 we show the SED of the source, 
including the values obtained at 1.25 mm and 2.22 fxm by Lefloch et al. (1996) and 2.65 mm 
by Ladd & Howe (1997). We have overplotted (following Ladd & Howe 1997) four simple gray 
body models, at Tdust = 10, 20, 30 and 40 K, of the form = B^{T) (1 - e"^) with 
T oc and the optical depth is normalized to the 2.65 mm flux. The simple models do at 
reasonable job at the IRAS and millimeter wavelengths for T^ust = 20 - 30 K, but are unable 
to fit the shorter wavelengths. 

It is possible to build more sophisticated models for the SED which include a density 
and temperature structure for the cloud core or envelope which surrounds the Class object 
(Andre & Montmerle 1994), and which also take into account the dust opacity as a function 
of chemical composition and wavelength (see Appendix). In Figure 7 we display three models 
with different inner gas cloud core densities (6 x 10^, 8x10^ and 10^ cm"^), but otherwise 
identical in their input parameters. The model with an initial density of 8x10"^ cm^^ (solid 
line) fits the observations amazingly well and is consistent with the upper limits at 2.22 //m and 
9.855 yum, which are based in 'no detections'. The models also illustrate how well they fit the 
IRAS and sub/millimeter observations, and how sensitive they are to the presence of the ~ 20 
fim. silicate feature. This feature essentially disappears in the 6 x 10^ cm~^ model (dotted line) 
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and not surprisingly becomes deeper with an increasing density at 10^ cm ^ (dashed hne). 

The model SEDs in Figure 7 assume a dust opacity dominated by bare sihcates, hence 
the strong absorption features at 10 and 20 //m, and a dust temperature of 18 K. The models 
consider a power law density and temperature distributions (see Appendix), with a core inner 
radius of 0.065 AU, and outer radius of 0.1 pc. The total masses for the cloud envelope 
surrounding the IRAS source for the three models (increasing in density) are 13, 17 and 22 
Mq. Our best model (with Tdust = 18 K and Menv = 17 Mq) yields different values than those 
obtained by using a constant density distribution and the simple gray body models, i. e. Tdust 
= 20 K and Menv = 10 M© (Ladd & Howe 1997). 



4. Conclusions 

We have analyzed some of the physical characteristics of Cep E embedded outflow using 
ISOCAM images in the v=0-0 S(5) 6.91 /xm and S(3) 9.66 fim molecular hydrogen lines. We 
find that the morphology of Cep E outflow of the ground vibrational H2 lines is similar to 
that of the near infrared emission in the v=l-0 2.12 ^m line. At these wavelengths, and at 
surface brightness of 12 - 46 /iJy/arcsec^, we do not detect the second H2 outflow almost 
perpendicular to the main 2.12 fj,m flow, nor we found traces of H2 along the second ^^CO 
J = 2 — 1 outflow at ~ 52° angle. 

We detect at 6.91 yum the likely source of the main H2 and CO outflows, IRAS 23011+6126, 
and show that the source is well detected in all IRAS bands using HiRes images. The source 
is not detected at 9.66 //m(nor at 9.54 fim), but we think this agrees with the interstellar 
extinction curve which has a minimum at ~ 7 //m, but rises a ~ 9.7 /im due to the strong 
silicate absorption feature, further enhanced in this case by a cocoon surrounding the Class 
object as the model of the SED seems to indicate. 
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The g^gj ratio is uniform and near unity across the outflow, a fact which is difficult 
to explain with simple plane parallel shock models (Eisloffel et al 1996). A constant S(5) to 
S(3) ratio also indicates that the extinction, with the exception around the IRAS 23011+6126 
source, is not affecting the H2 emission of the outflow lobes at these wavelengths. 

Assuming that the main source of opacity around IRAS 23011+6126 is due to bare 
sihcates, then our best model of SED for the envelope surrounding the Class source yields a 
total mass of 17 and a dust temperature of 18 K, 

We thank Ken Ganga for his help and insight on the ISOCAM calibration and the data 
reduction procedures. Our gratitude goes also to Jochen Eisloffel for sharing with us the 
analysis of Cep E with ISO. We thank N. King for obtaining the imaging data, and the referee 
for helpful comments and a careful reading of the manuscript. 
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Appendix 

The envelope around Cep E source is modeled as a series of spherically symmetric dust 
shells, where temperature and density are assumed to vary according to radial power laws 
whose exponents are free parameters. Once the external radius of the envelope is fixed, the 
inner envelope boundary is also determined and is equal to the radius where dust attains 
its sublimation temperature (~1500 K). As our aim was not only to provide a reasonable 
model for the global spectral energy distribution of Cep E source, but also to verify if the non 
detection at 9.67 jjm could be due to silicate absorption, we then decided to adopt the dust 
opacities tabulated by Ossenkopf & Henning (1994) for an MRN ( Mathis, Rumpl & Nordsieck 
1977) silicate dust with variable volumes of ice mantles, instead of the usual assumption (e.g. 
Hildebrand 1983) of a dust opacity simply expressed as a power law of frequency. 

Dust emission is computed from each shell as: 

F,,i = ^,,B{u,Ti)piVi (1) 

where k,, is the dust mass opacity, B{v,T) is the Planck function, pi is the density and 
Vj is the volume of i^^ shell. First, for each frequency the optical depth is computed inward 
as r = J2i{i^uPidfi) until either r = 1 is eventually reached for a certain shell, or the inner 
envelope radius is reached; having identified the inner shell which contributes to observed 
radiation, the emitted flux is computed outward according to Eq. |l] until the external radius is 
reached. The flux from each shell is extincted by the intervening shells toward the observer, 
and flnally the contributions from all shells are summed at each frequency. 
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Figure Captions 

Fig. 1. — (a) Grayscale image of the Cep E outflow in the v=0-0 S(3) 9.66 /im H2 hne (with 
continuum) in a 2' field of view, (b) The same grayscale image (continuum subtracted) with 
contours beginning at 0.4 mJy/arcsec^, and increasing by a factor 2^/^. 

Fig. 2. — (a) Grayscale image of the Cep E outfiow in the v=0-0 S(5) 6.91 //m H2 line (with 
continuum) in a 2' field of view, (b) The same grayscale image (continuum subtracted) with 
contours beginning at 0.2 mJy/arcsec^, and increasing by a factor 2^/^. 

Fig. 3. — Grayscale image of the v=0-0 S(5) 6.91 /xm line (with continuum) with a linear 
contour map of the v=l-0 2.12 jim emission. 

Fig. 4. — IRAS HiRes images at (from top left to right) 12, 25, 60 and 100 /xmcentered on the 
IRAS 23011+6126 source. The field is one degree, with the usual orientation, and the circle 
surrounding the IRAS source has a 5 ' radius. 

Fig. 5.— The v=0-0 S(3) 9.66 fj,m to v=0-0 S(5) 6.91 fj,m ratio across the Cep E outfiow. 
The ratio is essentially constant and near unity, except around the IRAS 23011+6126 source 
which is no detected at 9.66 /xm, and at the edges due to an artifact. 

Fig. 6. — The spectral energy distribution of IRAS 23011+6126 source (following Ladd & 
Howe 1997). The triangle corresponds to our 6.855 /im measurement, and the upper limits 
are for the 9.535 iim and 2.22 iim integrated fiuxes. The IRAS fiuxes were measured from the 
HiRes images. Overplotted for comparison four simple gray body models with T^ust = 10, 20, 
30 and 40 K, normalized at the 2.65 mm fiux 

Fig. 7. — More detailed models of the SED for IRAS 23011+6126 source using the same data 
as in Figure 6. The models assume a dust opacity dominated by dust grains made of bare 
sihcates, hence the strong absorption features at 10 and 20 //m (see the appendix for details). 
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The input parameters of the models are identical except for the initial inner core gas density 
of 6 X 10^ cm~^ (dotted line), 8 x 10^ cm~^ (sohd line) and 10^ cm~^ (dashed line). 
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